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We study dynamic ground state properties in the crossover from weak (BCS) to strong coupling 
(BEC) superfluidity. Our approach is based on the attractive Hubbard model which is analyzed 
by the dynamical mean field theory (DMFT) combined with the numerical renormalization group 
(NRG). We present an extension of the NRG method for effective impurity models to self-consistent 
calculations with superconducting symmetry breaking. In the one particle spectra we show quan- 
titatively how the Bogoliubov quasiparticles at weak coupling become suppressed at intermediate 
coupling. We also present results for the spin and charge gap. The extension of the NRG method to 
self-consistent superconducting solutions opens the possibility to study a range of other important 
applications. 
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Introduction - There has been a resurgence of inter- 
est into the nature of the crossover from the weak cou- 
pUng Bardeen, Cooper, and Schrieffer^ (BCS) superflu- 
idity to the Bose Einstein condensation (BEC) of pre- 
formed pairs, due to its recent experimental reahzation 
for ultracold atoms in an optical trap^'^'^'^. In the ul- 
tracold atom experiments, the interactions between the 
fermionic atoms in the trap can be tuned by a Fesh- 
bach resonance permitting an unprecedented control of 
the parameters, so that theoretical predictions may now 
be put to an experimental test. Above the transition 
temperature the two limiting cases, BCS superconduc- 
tivity and BEC, correspond to quite different physical 
situations. The weak coupling BCS theory describes a 
fermionic system with Cooper instability, whereas the 
strong coupling BEC limit is a system of strongly bound 
fermion pairs, which obey Bose statistics. The theoret- 
ical understanding which has been developed over the 
years is that the properties, such as the order parameter 
Asc and the transition temperature Tc to the superfluid 
state, are connected by a smooth crossover, and approx- 
imate interpolation schemes between these limits have 
been devised^i^i^. 

There is experimental evidence that the BCS-BEC 
crossover has also relevance for strong coupling and high 
temperature superconductors. These superconductors 
display properties, it has been claime d^°i^^i^^ , that can 
be understood in terms of local pairs, preformed above 
the transition temperature Tc, in contrast to the BCS 
picture, where the pairs no longer exist above Tc. 

As more experimental techniques, such as pairing gap 
spectroscopyi^ii^, are being developed to probe this 
crossover, more detailed theoretical predictions are re- 
quired. In particular the dynamic response functions 
have received little theoretical attention, and it is the pre- 
dictions for these quantities through the crossover that 
will be the focus of the present paper. 

Here, we use the attractive Hubbard modeU^ to study 
the crossover. In order to investigate the full crossover 
regime from weak to strong coupling a reliable non- 



perturbative method is needed. We employ the dynami- 
cal mean field theory (DMFT), together with the numer- 
ical renormalization group (NRG) to treat the effective 
impurity problem. DMFT studies, which are appropri- 
ate for the three and higher dimensional case, for the 
attractive Hubbard model with other methods for the 
effective impurity problem have been carried out in the 
normal phase^^^, and more recently in the symmetry- 
broken phaseiiiiiii^. But spectral properties have not 
been calculated and analyzed in detail there. There have 
been calculations of the dynamic response functions for 
the two dimensional model, for example with quantum 
Monte Carloi^ or recent work using the cellular DMFT—. 
Our results based on the DMFT are not expected to be 
directly applicable to this low dimensional case, but a 
comparison can be of interest to see if there are common 
trends. 

Our method to study the effective impurity problem, 
the NRG, is a non-perturbative method originally de- 
vised by Wilson?J to analyze the many body problem in 
the Kondo model and Anderson impurity model (AIM). 
In this approach the high energy degrees of freedom are 
eliminated numerically to find a an effective low energy 
theory, and its early application contributed vitally to the 
quantitative understanding of the Kondo probleirii^^. The 
method was developed further in the 1990s, when it was 
extended to the calculation of spectral functions^^i^i. Af- 
ter the advent of the DMFT— it was realized then that it 
could be used for self-consistent solutions of the effective 
impurity problem^^. The calculation of spectral func- 
tions was improved further by approaches based on the 
density matrix^l and the complete basis set proposed by 
Anders and Schille r^^i^^i^° . Our approach uses this latest 
full density matrix scheme to calculate reliable spectral 
functions. The NRG method and its applications are 
comprehensively reviewed by Bulla et ali^. 

In the original theory and in many of the applications 
one deals with a metallic system, where the NRG inher- 
ent logarithmic discretization and the descent to low en- 
ergies gives the accurate low energy properties. However, 
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the method has also given reUable results in situations 
with a gap in the spectrum. Examples include the Mott 
transitions^ or antiferromagnetic solution o^^i^"* in the re- 
pulsive Hubbard model in the DMFT-NRG framework. 
There is also an extensive amount of work for situations 
where the bath of the impurity model is a mean field BCS 
superconductor with given energy gap 35i36i37i38, 39,40,41 _ 

The NRG parameters for the bath in this case are fixed 
and can be determined by a straightforward generaliza- 
tion of the original NRG approach. Here, we present 
the extension of the method to calculations where the 
parameters describing the superconducting bath have to 
be determined self-consistently, which is a requirement of 
the DMFT approach, where the impurity is an effective 
one and plays an auxihary role. This opens the possi- 
bility to study a number of interesting non-perturbative 
problems involving superconductivity in a well controlled 
framework. 

Formalism - Our study is based on the attractive 
Hubbard modeU^, 



H = -Y^ iUjcl„Cj,a + h.c.) - ^ ^ - [/ ^ 



(1) 



with the chemical potential ^, the interaction strength 
U > and the hopping parameters Uj. cj^ creates a 

fermion at site i with spin a, and rii^a = cl^Ci^a-- For 
the DMFT calculation we employ the Anderson impu- 
rity model in a superconducting medium as the effective 
impurity model. 



HAnd = gjmp + Skcj ^Ck^a + Vkicj^^do- + h.C.) 



where iJimp = Y.a^d.nd^a - Und^fUd^ and is the 
fermionic operator at the impurity site. Efc, Vk and 
Afc are parameters of the medium. The non-interacting 
Green's function matrix at T = has the form, 



K_{lo) is the generalized matrix hybridization for the 
medium, with diagonal part 



and offdiagonal part. 



AD 



(4) 



(5) 



In the DMFT with superconducting symmetry breaking 
the effective Weiss field is a 2 x 2 matrix G_Q^{t)^. The 
DMFT self-consistency equation in this case also is a ma- 
trix equation. 



(6) 



with the fc-independent self-energy S(a;) and the local 
lattice Green's function G(w). We identify as usual 
Go{uj) — G_q{oj). Having calculated the local Green's 
function G and the self-energy the self-consistency 

equation |[6]) determines the new Weiss field and medium 
as input for the effective impurity problem. Once Kii{uj) 
and i^2i(w) are given by ([6]), the problem is to calculate 
the effective impurity model parameters and map ^ to 
the linear chain Hamiltonian, 

i^And = Himp + '^finifLafn+Ua +h.C.) (7) 

(7.71 

+ Y.^nfljn,a - J] A„ + h.C.) , 

(T,n n 

to which the iterative diagonaHzation of the NRG can 
be applied. e„, and A„ are the parameters of the 
linear chain model and fX ^ the operators for the sites. 
The details of how this can be achieved will be published 
elsewhere^. 

Results - Static and integrated quantities like the or- 
der parameter, the average pair density {n-^ni) or super- 
fiuid density Ds have been discussed in other works^^^. 
We will present our DMFT-NRG results for these quan- 
tities in a separate publication^. Here we focus on dy- 
namic response functions, which have received little at- 
tention so far. 

For numerical calculations within the DMFT-NRG ap- 
proach we take the semi-elliptical form of the Bethe lat- 
tice for the non-interacting density of states po(e) = 
2V£>2 - e'^liiD'^, where 2D is the band width with D = 
2t for the Hubbard model, t = 1 sets the energy scale 
in the following. All the results presented here are for 
T = and the generic case of quarter filling, n — 1/2. 

The spectral gap and the quasiparticle excitations 
can be analyzed from the one particle Green's function 
Gk{to), which is the diagonal element of the inverse of 
Gkiu^)-' 



(8) 



(3) In the BCS limit the excitations can be described by 



1 — with 



where ul = (1 + + Un/2)/E"^)/2, vl 
ik — £k — M- The mean field order parameter Asc.mf = 
U {cq,^Cq^i)m_f and the chemical potential /i are deter- 
mined from the gap and number equation, respectively. 
In this approximation the two bands of quasiparticle exci- 
tations are given by ±El. = ±yji^k + Un/2Y -\- A^^ ^p, 

with weights uj, for the positive and Vk for the negative 
excitations with infinite lifetime. In the real interacting 
system this is obviously not the case and our calculation 
of the dynamic quantities allows one to study this. 

The standard way to analyze the real quasiparti- 
cle excitations Ek is to use the many-body definition. 
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ReGfc(a; = Ek)~^ = 0. When the self-energy functions 
are calculated from DMFT this can be solved for a given 
Ek as an implicit equation. Due to the symmetries of the 
self-energy for a solution lu = also to — —Ek is a so- 
lution of this equation. In order to extract quasiparticle 
parameters in a simplified form, we expand around these 
solutions.— We define 



z{Ek)-^ = 



1 dReGkiuJ = Ek)-^ 



(9) 



2Ek duj 
In the vicinity oiuj = Ek the Green's function then reads 



u^Ek 



z{Ek)-Hu;-Ek)+iWiEk)' 



(10) 



where we have set uj — Ek in the numerator, v?{Ek) = 
[Ek + Cfe ~ ^22{Ek))/2Ek, and approximated the imagi- 
nary part by 



W{Ek) = 



ImGkji^ = Ek)-^ 
2Ek 



(11) 



independent of w. The main contribution for the spectral 
density near uj = Ek is then 



Pfc(w) = u^{Ek)z{Ek) 



W{Ek)/Ti 



{uj - Ekf + WiEkY 



(12) 



where W{Ek) = z{Ek)WiEk). If W{Ek) is small 
then this is a well-defined Lorentz quasiparticle peak 
with width W{Ek) and spectral weight w+{Ek) — 
z{Ek)v?{Ek)- The same analysis can be done near 
Lo — —Ek- With these quantities we can analyze the 
question up to which interaction strength there are well- 
defined fermionic (gapped) quasiparticles, i.e. which are 
not too broad. 

In Fig. [T] we give a typical example of efc-resolved 
spectra /3fc(ijj) for U — 2, which is the critical interaction 
Uc for bound state formation in the two-body problem 
for the Bethe latticei^. 

We can see the two bands of quasiparticle excitations 
with the sharpest peaks in the region of minimal spectral 
gap. These spectra can be compared with the ones which 
were calculated by Garg et al.-i^ by iterated perturbation 
theory. There the quasiparticle excitation delta peaks are 
disconnected from the continuum, which is however an 
artifact of the approximation for the self-energy, whose 
imaginary part vanishes over too large a region in uj. The 
mean field results (red arrows) describe the form of the 
quasiparticle bands qualitatively well. Also the weight of 
the peaks in the full spectrum pfc(tj) is comparable with 
the height of the arrows of the mean field theory. Notice 
that the sharpest peaks of pk{oj) have the Lorentzian 
shape as in (fT2l) . For different Ek and also at larger U 
the peak form can be asymmetric and (fT2|) not such a 
good fit. 

The effect of the dynamic fluctuations on the excitation 
spectrum can be seen clearer in the inset of Fig. (TJ where 




FIG. 1: (Color online) The Efe-resolved spectral functions 
pk{i^) for U = 2. The arrows show the delta-function peaks of 
the mean field result Pfc'^^(aj), where the height of the arrow 
indicates the weight of the peak. Inset: Mean field bands 
(dashed) and real quasiparticle band Ek (full line). 



we compare the mean-fleld bands E^ with the real quasi- 
particle band Ek- We can see a substantial reduction of 
the excitations energy for a certain bare energy Ek- For 
instance, the minimal spectral gap Acx — niin(£'fc) is re- 
duced by a factor of about 2.3 in this case. For different 
values of the interaction in the whole crossover regime 
this reduction effect can be seen in Fig. (2 where we 
plot the value of the order parameter computed by mean 
fleld (MF) theory, Asc,mf, and DMFT, Asc.dmft, and 
the spectral excitation gap Acx- Note that only at weak 
coupling the minimal spectral gap is strictly given by 
Acx, since at intermediate to stronger coupling there are 
broader peaks and one flnds flnite weight for \uj\ < Acx 
implying a smaller excitation gap.— 




FIG. 2: (Color online) Comparison of the order parameter 
Asc (mean field theory and DMFT) with spectral g for 
a range of interactions U- 

The order parameter Asc,dmft is reduced from its mean 
field value Asc,mf upon including fiuctuations for all cou- 
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pling strengths4^ Furthermore, we can see that the gap 
in the single particle spectrum Aox, as deduced from the 
one particle Green's function, is in turn smaller than 
Asc.DMFT for intermediate to strong coupling. At weak 
coupling these quantities tend to the same value suggest- 
ing an excitation spectrum of a renormalized mean field 
form. 

In the weak coupling mean field theory the elementary 
excitations are gapped fermions. However, for strong 
coupHng we deal with a bosonic system, and due to the 
large gap of order U fermionic excitations are suppressed. 
We now analyze at which point in the intermediate cou- 
pling regime one tends to lose the sharp quasiparticle ex- 
citations. In order to study this question we look at the 
width of the excitations defined at the smallest spectral 
excitations W — W{Ek — Aex). If this is small we have 
a strongly peaked well-defined quasiparticle excitation as 
observed in Fig. [TJ whereas for larger width the excita- 
tion decays more rapidly. For different local interaction 
strengths this is illustrated in Fig. [31 




FIG. 3: The quasiparticle width W at the excitation gap 
energy for a range of interactions U. Inset: The ratio W/Acx- 

We find that whilst W is very small for weak coupHng 
it increases rapidly in the intermediate coupling regime, 
U = 2-4. The ratio W/Ac^, which is plotted in the 
inset of Fig. [31 also increases from a small value at weak 
coupling to larger values at intermediate coupling. This 
shows that sharp (long-lived) quasiparticles excitations 
are suppressed then. 

We now turn to spin and charge response functions. 
Within our framework we can calculate the local dy- 
namic spin susceptibility Xs(w) = i^/N)J2qXs{'^,Q) and 
charge susceptibility Xc(w)- The temperature depen- 
dence of Xs in the static limit has been discussed by 
Keller et al.~ and the existence of a spin gap was demon- 
strated. In the following Fig. [H we show the imaginary 
part of Xsi^) (top) and Xc(w) (bottom) for a number of 
values of U. 

The spin response is suppressed when the local attraction 
U is increased. The spin gap^ Agp grows with U as can 
be seen in the inset of Fig. \E[ and it is directly related 
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FIG. 4: (Color online) The imaginary part of the local dy- 
namic spin susceptibility Xs (i^) (top) and of the charge sus- 
ceptibility Xc('^) (bottom) for different U . 



to the gap in the one particle excitation spectrum. This 
is expected since spin excitations are accompanied by 
breaking of a fermionic pair with opposite spin. For large 
U the binding energy is proportional to U. The charge 
susceptibility shows a different behavior. The broad peak 
for weak coupling becomes sharper in the strong coupling 
limit and moves to lower energy. A similar trend is seen 
in the results for the charge susceptibility for q = in the 
two dimensional model^S. A peak related to the massless 
Goldstone mode ~ is not present in our DMFT 
calculation for the local charge susceptibility. The vis- 
ible charge gap^ Ac, which is always smaller than the 
spin gap, first increases with U but at strong coupling 
decreases again (see Fig. [5]). The maximum occurs at a 
value, which is a bit larger than Uc = 2 for the two-body 
bound state problem. 

One way of understanding this behavior is to consider 
the mapping of the attractive Hubbard model to the half 
filled repulsive model in a magnetic field. At strong cou- 
pling this can be mapped to a Heisenberg model with 
spin coupling J = 4t'^/U. Charge excitations in the at- 
tractive model correspond to spin excitations in the re- 
pulsive model, and the characteristic scale for the latter 
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FIG. 5: (Color online) The charge gap Ac for a range of 
interactions U . The dashed line behaves ~ At^ /U. Inset: The 
spin gap as a function of U. 

is J. The dashed hne in Fig. [5] proportional to At'^/U 
with an adjusted prefactor reflects this dependence. 

Perspective - Let us comment on the range of apphca- 
tions of the extended DMFT-NRG method. Within this 
approach low temperature static and dynamic response 
functions can be calculated. In the experimental situa- 
tion for the BCS-BEC crossover, the trapped fermionic 
gases constitute an inhomogeneous system. With a 
larger numerical effort these situations can be studied 
with real space DMFT extensions as done for the Mott 
transition^. In the experiments it is difficult to cool the 
fermions down to very low temperatures compared to 
the Fermi energy. At flnite temperature, but still in the 
symmetry-broken phase, the excitation gap is expected 
to be smaller compared to our T — case and the width 
of the quasiparticle excitation largeriS.. A calculation at 
flnite temperature would be desirable to see how well 
the general features described here are still visible in the 
whole crossover regime. 

Moreover, with the given method other models acces- 
sible by DMFT can be investigated. This includes, for 
instance, systems with a local electron phonon coupling, 
which are relevant to understand the fulleride supercon- 
ductors. Also periodic Anderson type models and occur- 



ring superconducting phases are accessible. 

Apart from these approximation to lattice models 
with a local self-energy, our method could also be 
used to study open problems for magnetic impurities 
in superconductors^^. According to Anderson's theo- 
rem and the work of Abrikosov and Gor'koA^Sii^ a fl- 
nite concentration of non-magnetic impurities does not 
alter an s-wave superconducting state, but magnetic im- 
purities suppress the order substantially. So far there 
are, however, few reliable quantitative microscopic stud- 
ies of the local effect of a single impurity (magnetic or 
non-magnetic) on a superconductor. Since this is a non- 
perturbative problem our self-consistent NRG calculation 
could be an excellent method to tackle this problem. An 
attempt to investigate this based on the NRG method 
without self-consistency was made by Sakai et al.— . 

Conclusions - We have applied an extended DMFT- 
NRG method to calculate self-consistent solutions with 
superconducting symmetry breaking in the attractive 
Hubbard model. This method can access static and dy- 
namic quantities for any coupling strength and doping. 
We have focused on the zero temperature one particle 
spectrum and the local dynamic spin and charge sus- 
ceptibilities for quarter fllling. It was shown that the 
excitation energies are reduced substantially from their 
mean fleld results due to fluctuations. At weak coupling 
there are well-deflned fermionic quasiparticle excitations, 
but as displayed in Fig. [3] the width of these excita- 
tion increases at intermediate coupling, such that sharp 
quasiparticle excitations are suppressed. In addition, the 
different behavior of the excitation gaps in the dynamic 
charge and spin susceptibilities has been analyzed quan- 
titatively. We also pointed out that this work is not only 
relevant for the attractive Hubbard model, but the ap- 
proach can be used to study superconductivity in other 
models, as well as for a microscopic description of the 
effect of impurities on superconductors. 
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